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ABSTRACT 



Context. Numerous spectroscopic observations provide compelling evidence for non-canonical processes that modify the surface 
abundances of low- and intermediate-mass stars beyond the predictions of standard stellar theory. 

Aims. We study the effects of thermohaline instability and rotation-induced mixing in the 1-4 M ;, range at solar metallicity. 
Methods. We present evolutionary models by considering both thermohaline and rotation-induced mixing in stellar interior. We 
discuss the effects of these processes on the chemical properties of stars from the zero age main sequence up to the end of the second 
dredge-up on the early-AGB for intermediate-mass stars and up to the AGB tip for low-mass stars. Model predictions are compared 
to observational data for hthium, '^c/'^c, [N/C], [Na/Fe], i^O/'^O, and i''0/"*0 in Galactic open clusters and in field stars with 
well-defined evolutionary status, as well as in planetary nebulae. 

Results. Thermohaline mixing simultaneously accounts for the observed behaviour of '^C/'-'C, [N/C], and lithium in low-mass stars 
that are more luminous than the RGB bump, and its efficiency is increasing with decreasing initial stellar mass. On the TP-AGB, 
thermohaline mixing leads to lithium production, although the ^Li yields remain negative. Although the ^^He stellar yields are much 
reduced thanks to this process, we find that solar-metallicity, low-mass stars remain net "'He producers. Rotation-induced mixing is 
found to change the stellar structure so that in the mass range between ~ 1.5 and 2.2 M,j, the thermohaline instability occurs earlier 
on the red giant branch than in non-rotating models. Finally rotation accounts for the observed star- to-star abundance variations at a 
given evolutionary status, and is necessary to explain the features of CN-processed material in intermediate-mass stars. 
Conclusions. Overall, the present models account for the observational constraints very well over the whole mass range presently 
investigated. 
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1. Introduction 

At all stages of their evolution, low- and intermediate-mass 
stars (LIMS) exhibit the signatures of complex physical pro- 
cesses that require challenging modelling beyond canonical (or 
standard) stellar theorjQ. The combined effects of rotation- 
induced mixing, atomic diffusion, and internal gravity waves, 
have been extensively studied and were shown to account self- 
consistently for observational patterns of light elements such 
as lithiu m and beryllium in rnain sequence and subgiant L IMS 
(see e.g. ICharbonnel & TalonI l200"8l Ismilianic et alj|2009l and 
refere nces therein). During the first dredge-up (IDUP; llbeni 
Il967h . the stellar surface chemical composition is further mod- 
ified when the deepening convective envelope mixes the ex- 
ternal layers with hydrogen-processed material. Dilution then 
changes the surface abundances of helium-3, lithium, beryl- 
lium, boron, carbon, nitrogen, and eventually sodium, with mod- 
ification amplitudes that strongly depend on the initial stel- 
lar mass and metallicity (e. g. Sweig ait et alj|1989l : ICharbonnell 
11994 1: 'Boothroyd & Sackmann 1999). After the completion of 
the first dredge-up, both standard and rec ent rotating models 
jChaname et alj 120051: iPalacios et"al] |2006|) predict no further 



variations in the surface abundance patterns until the stars start 
climbing the asymptotic giant branch (AGB). 

Numerous observations provide, however, compelling 
evidence of a non-canonical mixing process that occurs when 
low-mass stars reach the so-called bump in the luminos- 
ity function on the red giant branch (RGB). At that phase, 
indeed, the surface carbon isotopic ratio drops, together 
with the abundances of lithi um and carb o n, whi le that of 



nitro 



en 



increase s slightly jGilrovl 1 19891: Gilrov & Brown 
11991; 'Luck' '1994; 'Chai-bo nne]|ll994t Oiarbonnel et al. 199"! 
Charbonnel & Do NascimenT oT Il998t iGratton et al. 200C 



Tautvaisiene et al. 2000; TaiTtvaisie ne et al.l 120051; ISmith et al 



2002 - Shetronel i2003c iPilachowski et al.1 120031 



IqdS LSpit eet aLl bTM 



Geisler et al 

■■ .- .Recio-Blanco & de Lavernvl 120071 ; 
Smilianic et al. '2009^. Based on lithium observations, 
iCharbonnel & Balachandran (.2000) proposed that intermediate- 
mass stars suffer from a similar extra-mixing episode when they 
reach the equivalent of the bump on the early-AGB phase, after 
helium exhaustion. 

Thermohaline mixing has recently been identified as the 
mechanism that govern s the photospheric compo sition of 
low-mass brigh{3 giants (ICharbonnel & ZahnI l2007bl hereafter 
CZ07). In such stars, this double diffusive instability is in- 



' By canonical we refer to the modelling of non-rotating, non- 
magnetic stars, in which convection is the only mechanism that drives 
mixing in stellar interiors. 



^ "Bright" refers here to stars that are more evolved than the RGB 
bump. 
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duced by the molecular weight inversion created by the 
^He(^He,2p)'*He reaction in the external wing of the hydrogen- 
burning shell. Indeed this peculiar reaction converts two parti- 
cles into three and thus d ecreases the m ean molecular weight, 
as already pointed out by lUMchl (1197 ll) although in a different 
stellar context. The thermohaline instability is expected to set in 
after the first dredge-up when the star reaches the RGB lumi- 
nosity bump. In terms of stellar structure, the RGB bump corre- 
sponds to the moment when the hydrogen-burning shell (here- 
after HBS) encounters the chemical discontinuity created inside 
the star by the convective envelope at its maximum extent during 
the first dredge-up. When the source shell (which provides the 
stellar luminosity on the RGB) reaches the border of the H-rich 
previously mixed zone, the corresponding decrease in molec- 
ular weight of the H-burning layers induces a drop in the to- 
tal stellar luminosity, thereby creating a bump in the luminosity 
func tion since stars spend a relatively long tim e at this location 
(i.e..lFusi Pecci et alj|l990l:ICharbonnellll994tlCharbonnel et al.l 



T999l:lBania et allll997il20Q2l:lTosilll998l:lDearborn et al.lll996l: 
Charbonnelll2002l:lRomano et al.ll2003ir ^ 



11998 ). Afterwards H-burning occurs in a region of uniform com- 
position, allowing for the molecular weight inversion due to ^He 
burning to show up and thus enabling the thermohaline instabil- 
ity to set in. 



Actually it was lEggleton et alj ( l2006l) who first drew atten- 
tion to the destabilizing role of the mu-inversion due to ^He- 
burning in a red giant. Their claim was based on a 3D sim- 
ulation aimed at studying the core helium flash of low-mass 
red giants with the hydrodynamic code "Djehuty" and per- 
formed with a mu-profile dr awn from a ID evolutionary se- 
quence ( iDearborn et al.ll2006l) . which fortuitously demonstrated 
what the authors called an "unexpected mixing" ; they ascribed 
it to the wel l-known Rayleigh-Taylor instability dEggleton et alJ 
120061 l2007h . However CZ07 pointed out that in a star, as the 
inverse mu-gradient is gradually building up, the first insta- 
bility to occur and to modify the mu-profile is the double- 
diffusive instability known in the Uterature u nder the generic 
name of thermohaline instability (ISternlll960h . This important 
precision on the actual nature of th e physical proc e ss ope rat- 
ing in the star was acknowledged by Eggleton et"al] ( 20081 see 



also Denissenkov & Pinsonneaultll20()8l and Cantiello & LangeJ 
120 id " it is not just a question of semantics, since these two 
instabilities - one dynamical and the other double diffusive 
- proceed on much different timescales. It is important to 
note that thermohaline instability has long been known to de- 
velop in other stellar situations whenever inverted molecular- 
weight gradients are built. This is the case for instance when 
helium- or carbon-rich material is de posited at the surface of 
a star in a mass-tran sferring binary ( IStothers & SimonI 119691 : 
IStancUffe et ani2007h . when a star a ccretes heavy elements dur- 
ing planet formation (IVauclairl2004l) . or after the ignitio n of ^He 
burning in a degenerate shell (i.e., core helium flash; see lThoniasI 
Il967lll970l) . It also develops in stars where radiative levitation 
lead s to the accurnulatio n of heavy elements in outer stellar lay- 
ers dTheado et al.ll2009l) . Since the term "thermohaline" is long 
established in the literature, we see no reason to replace it by the 
expression "5/i -mixing" as proposed bvFEggl eton et al. (20081). 

CZ07 showed how nicely models including thermohaline 
mixing as des cribed by the theoretical prescription advocated by 
lUMchl d 19721 see §2.2) do account for the carbon isotopic ratio, 
as well as for the lithium, carbon, and nitroge n abundances in 
low-m etallicity, low-mass bright giants (see also lStancliffe et aP 
l2009l) and simultaneously reduce significantly the stellar yields 
of ^He with respect to canonical models, as required by mea- 
surements of ^He/H in Galactic HII regions dBalseret al.lll994 



Note that the connection between the so-called "helium- 
3 problem" and the behaviour of the carbon isotopes in 
RGB stars was long established. Indeed although it had not 
yet been identified, the mechanism responsible for the low 
'^C/'-'C values in bright giants was expected to lead to the 
destruction of ^He by a large factor i n the bulk o f the stel- 
lar e nvelope as initially suggested b y Roodetal] (1984 see 
also iHoganI 119951 ICharbonnell ll995UWasserburg et all 11995 



-?ll995t 

We iss et al.lll996l;ISackmann & Boothrovdll999HEggleton et al.1 



2006'; 'Balser et al.' '2001). CZ07 results were confirmed by 
lEggleton et al.i (i2008i) although with a simplistic phenomenolog- 
ical procedure to est i mate the dif fusion coefficient. Additionally, 
ICantiello & Langeil d2008ll2010h reported that thermohaline in- 
stability can still occur during core helium-burning and beyond 
in stars that have kept a ^He re servo ir at these ad vanced phases. 
Finally IStancliffe et all d2009l) and ISiessI d2009h discussed the 
impact of this mechanism in low-metallicity intermediate-mass 
thermally pulsing AGB stars and in super-AGB stars respec- 
tively. 

The signatures of thermohaline mixing induced by ^He- 
burning have been observed in giants of both open and globular 
clusters, as well as in field stars including extremely metal-poor 
giants, and in stars belonging to extragalactic systems such as 
the Large Magellanic Cloud and the Sculptor galaxy (see refer- 
ences above). It appears thus to be a universal process that occurs 
independently of the stellar environment, although it may be in- 
hibited in so me rare RGB stars hosting st rong fossil magnetic 
fields. Indeed ICharbonnel & ZahnI d2007ah examined the effect 
of a magnetic field on the thermohaline instability, and con- 
cluded that in a large fraction of the descendants of Ap stars ther- 
mohaline mixing should not occur. The relative number of such 
stars with respect to non-magnetic objects that und ergo thermo- 
haline mixin g is very low (le ss than 5%, see i.e., IWolfjlI968l: 
NorthJ993MPower et al.l2007l) and consistent with the statistical 
constraint coming from observati ons of the carbon isotopic ratio 
in evo lved stars as estimated by ICharbonnel & Do Nascimentol 
(119981) . It also reconciles the measurements of ^He/H in Galactic 
HII regions with high values of ^H e observed in a couple of plan - 
etary nebulae (for more details see ICharbonnel & Zahnll2007ah . 

However, the impact of thermohaline mixing on the stellar 
chemical properties appears to depend on the initial stellar mass 
and metallicity, as suggested by the observations of the carbon 
isotopic ratio (which is the most reliable signature of this mech- 
anism) in giant stars over broad ranges in mass and metallicity 
(see references above). Until now, the theoretical background for 
these correlations is not firmly established (e.g., the mass range 
where thermohaline convection actually modifies the stellar sur- 
face composition), and the influence of other quantities such as 
the rotation history has never been investigated. 

This is the aim of this series of papers where we investigate 
the occurrence and the impact of thermohaline mixing in stars of 
various initial masses and metallicities with non-canonical mod- 
els, i.e., taking into account self-consistently this mechanism to- 
gether with rotation-induced mixing. Here we focus on the solar 
metallicity case, with the main observational constraints com- 
ing from abundance determinations in evolved stars with well 
defined evolutionary status up to the early-AGB and belonging 
to open clusters with turnoff masses between ~ 1 and 4 
or to the field. We also use as observational constraints lithium 
abundance determinations in low-mass oxygen-rich AGB vari- 
ables and carbon isotopic ratios in planetary nebulae. In § 2 we 
present the input physics of our models. In § 3 we discuss the 
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predictions of our models computed with thermohaline mixing 
only, and then of our models that also include rotation-induced 
mixing. Theoretical predictions are compared to observations in 
§ 4 before we conclude in § 5. 



2. Input physics for the stellar models 

We present models computed with the c ode STARE VOL 
dSiess et^ 1 120001: iPalacios et al.ll2003L l2006l) at solar metallic- 
ity (with Asplund et al. 2005 chemical c omposition except for 
Ne for which we use the value derived bv lCunha et al.ll2006l) for 
a range of initial stellar masses between 1 and 4 Mq. In order to 
quantify precisely the impact of each transport process at the var- 
ious evolutionary phases, we will show predictions for models 
computed with the following assumptions: (1) Standard mod- 
els (no mixing mechanism other than convection); (2) Models 
including thermohaline mixing only (rotation velocity V=0); 
(3) Models including thermohaline mixing and rotation-induced 
processes for different initial rotation velocities. All the mod- 
els are computed up to the end of the second dredge-up on the 
early-AGB. Some of the low-mass models are followed along 
the TP-AGB up to the AGB tip. 



2.1. Microphysics 



We use the OPAL opacity tables ( llglesias &Rogersl[T996h for 
T > 8000 K that acc ount for C and O enrichments, and the 
iFerguson et al.l (l2005h data at lower temperatures. In both cases 
tables consistent with the assumed initial composition have been 
generatecQ. We follow the evolution of 53 chemical species from 
to 3^C1 us ing th e nominal NACRE nuclear reaction rates 
jAngulo et al.l 1 19991) by default and those given in Appendix 
A otherwise. The equation of state is described in details in 
ISiess et al ] (I2000h and accounts for the non-ideal effects due 
to Coulomb interactions and pressure ionisation. The treatment 
of convection is based on the classical mixing length formal- 
ism with aMLT = 1-6, and no convective overshoot i s incl uded. 
The mass loss rates are computed with iReimersI (Il975h for- 
mula (with rjR = . 5) up to the early-AGB phase, and with 
IVassiliadis & Wood (Il993h prescription on the TP-AGB. 



2.2. Thermohaline mixing 

The thermohaline instability occurs in a stable stratification that 
satisfies the Ledoux criterion for convective instability: 



Vad-V- 



Vm>0, 



but where the molecular weight decreases with depth: 

V,:=^<0 
^ dlnP 



(1) 



(2) 



with the classical notations for V = ((9lnT/(9lnP), (p = 
((9lnp/(9lnjLi)p,T and 5 = — ((9lnp/5lnT)p,^, and Vad be- 
ing respectively the molecular weight gradient and the adiabatic 
gradient. 

For the turbulent diffusivity produced by the thermohaline 
instablity we use the prescription advocated by CZ07 based on 
[Ulrich (1972) arguments for the aspect ratio a (length/width) 
of the salt fingers as supported by laboratory experiments 



iKrishnamurtil (l2003h and including lKippenhahn et al.l (Il980l) ex- 
tended expression for the case of a non-perfect gas (including 
radiation pressure, degeneracy): 



Dt^QKl'I'l 77^-^ forV„<0, 



.5J (Vad-V) 
with K the thermal diffusivity. 



Ct = -TT^a^, 



(3) 



(4) 



and with a = 5 (Ulrich 1972) this coefficient is Ct=658. For con- 
sistency reasons we assume actually Ct=1000 as in CZ07. 

2.3. Rotation-induced mixing 

For the treatment of rotation-induced mixing we proceed as fol- 
lows. Solid-body rotation is assumed when the star arrives on the 
zero age main sequence (ZAMS). Typical initial (i.e., ZAMS) ro- 
tation velocities are chosen depending on the stellar mass based 
on ob served rotation distributions in young open clusters (Galgj 
Il993h . Surface braking by a magnetic torque is applied for stars 
with an effective temperature on the ZAMS lower than 6900 K 
th at have relatively a t hick c onv ective envelope as discusse d 
in lTalon & Charbonnell(ll998h and lCharbonnel & Talonl d 19991): 
the ad opted braking law follows the description of Kawaleil 
(119881) . From the ZAMS on the evolution of the internal angular 
momentum p rofile is acco unte d for with the complete formalism 
developed bv lZahnI (1 1 992h and lMaeder & Zahnl(ll998h that takes 
into account advection by meridional circulation and diffusion 
by shear turbulen ce (for a description of the implementatio n in 
STAREVOL. see iPalacios et al.ll2003l iPalacios et al]l2006l and 



^ http://opalopacity.llnl.gov 
and |http://webs. wichita.edu/physics/opacityl 



iDecressin et al .I2009h . The transport of chemicals resulting from 
meridional circulation and both horizontal and vertical turbu- 
lence is computed as a diffusive process throughout the evolu- 
tion. The complete treatment for the transport of angular mo- 
mentum and chemicals is applied up to the RGB tip or up to 
the second dredge-up for the stars with initial masses below or 
above 2.0 M© respectively. The convective envelope is supposed 
to rotate as a solid body (uniform angular velocity) throughout 
the evolution; we discuss the implications of this assumption in 
§ 3.2.2. The transport of angular momentum by internal grav- 
ity waves (which is efficient only in main sequence stars with 
effective temperatures on the ZAMS lower than 6500 K, see 
iTalon & Charbonnell2003l) . is neglected. 

In the present work the transport coefficients for chemi- 
cals associated to thermohaline and rotation-induced mixings are 
simply added in the diffusion equation and we do not consider 
the possible interactions between the two mechanisms, nor with 
magnetic diffusion. As a matter of fact under the present assump- 
tions the thermohaline diffusion coefficient is several orders of 
magnitude higher than the total diffusion c oefficient associated 
to rota tion (see §3). This is confirmed by ICantiello & LangeJ 
(l20Tol) who also show that magnetic diffusion in RGB stars 
is much less efficient than thermohaline mixing. However, we 
should keep an eye on future hydrodynamic calculations that are 
required to evaluate with confidence the possible interactions 
of thermohaline fingers with differential rotation and magnetic 
fields in red giants. 

3. Theoretical predictions 

We first consider the case of low-mass stars that ignite helium- 
burning by a flash at the tip of the RGB well above (in terms 
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Table 2. Surface values of '^c/'^c, [C/Fe], [N/Fe], [Na/Fe], N(Li) and N(Be) at the end of the first and second dredge-up (IDUP 
and 2DUP respectively) and at the RGB tip (RGB) for models computed under different assumptions: Standard (st, no thermohaline 
nor rotation-induced mixing), with thermohaline mixing only (th), and with both thermohaline and rotation-induced mixing (th-nrot). 



M 




^ZAMS 




i^c/'-<c 






[C/FeJ 






[N/Fe] 






N(Li) 






N(Be) 






[Na/FeJ 




(Mo) 




(faii.s-') 


IDUP 


RGB 


2DUP 


IDUP 


RGB 


2DUP 


IDUP 


RGB 


2DUP 


IDUP 


RGB 


2DUP 


IDUP 


RGB 


2DUP 


IDUP 


RGB 


2DUP 


1 


St 





28.8 


28.8 


25 


-0.05 


-0.05 


-0.07 


0.16 


0.16 


0.21 


1.1 


1.1 


0.85 


0.34 


0.34 


0.15 













th 





28.8 


8.1 


7.7 


-0.05 


-0.09 


-0.10 


0.16 


0.25 


0.28 


1.1 


-1.3 


-1.87 


0.33 


-1.65 


-2.1 











1.1 


th 





27.35 


8.9 


8.5 


-0.07 


-0.1 1 


-0.12 


0.22 


0.28 


0.33 


1.34 


-0.58 


-0.97 


0.27 


-1.33 


-1.63 











1.25 


St 





25.74 


25.78 




-0.10 


-0.10 




0.27 


0.27 




1.46 


1.46 




0.18 


0.18 














th 





25.6 


10.4 


9.9 


-0.10 


-0.12 


-0.14 


0.37 


0.31 


0.33 


1.45 


0.07 


-0.26 


0.18 


-0.97 


-1.23 













th+rot 


50 


23.69 


10 




-0.12 


-0.14 




0.39 


0.34 




-4.35 


-3.25 




-2.12 


-3.24 














th+rot 


80 


21.6 


9.6 


9.1 


-0.13 


-0.15 


-0.16 


0.31 


0.35 


0.36 


-6.2 


-3.2 


-3.9 


-3.7 


-4.8 


-5.3 













th+rot 


1 10 


18.6 


9.1 


8.5 


-0.14 


-0.16 


-0.18 


0.33 


0.37 


0.38 


-6.87 


-3.25 


-3.94 


-5.12 


-6.21 


-6.74 











1.3 


th 





24.8 


1 1 




-0.11 


-0.13 




0.39 


0.32 




1.46 


0.22 




0.15 


0.88 












1.4 


St 





24.1 


24.1 


23.3 


-0.13 


-0.13 


-0.14 


0.31 


0.31 


0.32 


1.46 


1.46 


1.37 


0.01 


0.01 


0.03 













th 





24 2 


12 14 




-0 13 


-0 14 




0.31 


0.33 




1 46 


0.46 




01 


-0 73 




Q 


Q 








Q 


23 2 


23 2 


2"t g 


-0 14 


-0 14 


-0 15 


33 


33 


34 


1 45 


1 45 


1 37 


06 


06 


^ 




Q 


Q 




th 





23.3 


15.3 




-0.14 


-0.15 




0.33 


0.34 




1.45 


1.03 




0.06 


-0.30 


_ 














1 10 


21 2 


12 8 


1 1 8 


-0.16 


-0.18 


-0.18 


36 


38 


39 


0.43 


-0 28 


-0 87 


-0 47 


-1 06 


-1.49 


0.02 


0.02 


02 


— — 




Q 


22 5 


22 5 


21 9 


-0 16 


-0 16 


-0 17 


36 


36 


37 


1 41 


1 41 


1 34 


-0 0"* 


-0 02 


-0 07 


Q 


Q 


Q 




th 





22.6 


18.6 


16.7 


-0.16 


-0.17 


-0.17 


0.36 


0.36 


0.37 


1.42 


1.24 


0.95 


-0.004 


-0.16 


-0.39 











1.8 


th 





22.2 


19.9 




-0.17 


-0.18 




0.37 


0.37 




1.41 


1.32 




-0.03 


-0.11 














th+rot 


1 10 


20.04 


15.24 




-0.18 


-0.19 




0.31 


0.31 




0.59 


0.19 




-0.42 


-0.75 




0.15 


0.15 




1 9 


th 





21 9 


20 2 


18 7 


-0.18 


-0 18 


-0 19 


38 


0.38 


39 


1.39 


1.31 


1.05 


-0 062 


-0 13 


-0 32 





Q 


Q 


2.0 


St 





21.8 


21.8 


21.3 


-0.19 


-0.19 


-0.19 


0.39 


0.39 


0.40 


1.38 


1.38 


1.31 


-0.07 


-0.07 


-0.11 


0.04 


0.04 


0.04 




th 





21.8 


20.8 


19.6 


-0.19 


-0.19 


-0.20 


0.40 


0.39 


0.40 


1.37 


1.31 


1.10 


-0.07 


-0.12 


-0.28 













th+rot 


110 


19.4 


16.6 


16.2 


-0.22 


-0.22 


-0.23 


0.47 


0.47 


0.47 


0.55 


0.33 


0.16 


-0.48 


-0.65 


-0.77 


0.18 


0.18 


0.18 




th+rot 


180 


17.7 


15.3 


15 


-0.23 


-0.23 


-0.24 


0.47 


0.48 


0.48 


-0.33 


-0.54 


-0.92 


-0.86 


-1.03 


-1.21 


0.18 


0.18 


0.18 




th+rot 


250 


14.8 


13.4 


13.2 


-0.22 


-0.23 


-0.23 


0.47 


0.48 


0.48 


-0.96 


-1.15 


-1.3 


-1.29 


-1.44 


-1.55 


0.22 


0.22 


0.22 


2.5 


St 





21.2 


21.2 


21.1 


-0.21 


-0.21 


-0.21 


0.45 


0.45 


0.45 


1.34 


1.34 


1.25 


-0.11 


-0.11 


-0.16 





0.16 


0.16 




th 





21.2 


21.4 


20.9 


-0.21 


-0.27 


-0.27 


0.45 


0.45 


0.45 


1.36 


1.36 


1.27 


-0.01 


-0.01 


-0.14 





0.16 


0.16 




th+rot 


300 


17.2 


14.7 


13 


-0.53 


-0.53 


-0.53 


0.55 


0.82 


0.86 


-4.99 


-4.99 


-5.12 


-3.77 


-3.77 


-3.91 





0.76 


0.85 


2.7 


th+rot 


110 


19.3 


19.3 


19.1 


-0.24 


-0.24 


-0.24 


0.52 


0.52 


0.53 


0.53 


0.53 


0.44 


-0.49 


-0.49 


-0.54 


0.28 


0.28 


0.28 




th+rot 


250 


16.4 


16.4 




-0.28 


-0.28 




0.56 


0.56 




-0.95 


-0.95 




-1.28 


-1.28 




0.32 


0.32 






th+rot 


300 


18 


15.1 




-0.53 


-0.53 




0.83 


0.84 




-4.9 


-4.9 




-3.79 


-3.79 




0.76 


0.76 




3.0 


th 





20.85 


20.84 


20.55 


-0.21 


-0.21 


-0.22 


0.47 


0.48 


0.48 


1.30 


1.30 


1.15 


-0.15 


-0.15 


-0.23 


0.031 


0.22 


0.23 


4.0 


St 





20.48 


20.48 


19.94 


-0.21 


-0.21 


-0.22 


0.49 


0.49 


0.50 


1.28 


1.28 


1.02 


-0.16 


-0.16 


-0.33 





0.26 


0.28 




th 





20.49 


20.49 


19.86 


-0.20 


-0.20 


-0.22 


0.49 


0.49 


0.51 


1.29 


1.29 


1.03 


-0.16 


-0.16 


-0.33 


0.25 


0.25 


0.28 




th+rot 


300 


14 


14 


13.82 


-0.31 


-0.31 


-0.32 


0.61 


0.61 


0.62 


-1.76 


-1.77 


-1.94 


-1.77 


-1.77 


-1.87 


0.4 


0.4 


0.42 



Table 1. Luminosities of the bump (Lb min and Lb max. see Fig[T]), 
of the evolution point when the thermohaline zone contacts the 
convective envelope (Lc th), and of the RGB tip (L^p), for the 
low-mass stars at various initial rotation velocities. 



M 


yzams 


Lb. min 


Lb, max 


Lc.th 


Ltip 


(Ms) 




(U,) 


(Lq) 


(La) 


(Lg) 


1.25 





37 


45 


94 


2821 




50 


25 


25 


64 


2798 




80 


15 


16 


72 


2807 




110 


16 


17 


72 


2798 


1.5 





50 


59 


784 


2903 




110 


27 


36 


101 


2869 


1.8 





74 


83 


1907 


2995 




110 


51 


57 


145 


2768 


2.0 





79 


87 


1908 


2994 




110 


64 


87 


250 


2872 




180 


74 


117 


256 


2670 




250 


76 


95 


232 


2416 







3.1. 1.25 M0 model with thermohaline mixing only 

Figure [T] presents the evolutionary track in the Hertzsprung- 
Russell diagram (HRD) of the 1.25 M© model computed with 
thermohaline mixing only (no rotation). Several points are se- 
lected along the track in order to discuss the evolution of some 
relevant stellar properties. Aj 25 corresponds to the turnoff, when 
the hydrogen mass fraction in the stellar core is below 10^**. 
Bi.25 is chosen at intermediate luminosity between the bump 
(which minimum and maximum luminosities, Lb min and Lb max. 
are also shown) and the moment when the thermohaline zone 
"contacts" the convective envelope (see below). C1.25 stands at 
the "contact" luminosity L^.th where surface abundances start 
changing due to thermohaline mixing. Di 25 and Ei 25 are close 
to and at the tip of the RGB (then the mass of the helium core is 
0.428 and 0.486 Mq respectively, for a total stellar mass of L14 
and L03 Mg). 



3.1 .1 . Main sequence and subgiant branch 



of luminosity) the bump. With the considered metallicity and in- 
put physics this corresponds to stars with initial masses below 
^2.2 Mq. We present detailed predictions for a 1.25 M0 model 
computed without and with rotation (but with thermohaline mix- 
ing in both cases) in § 3.1 and 3.2 respectively, and discuss the 
uncertainties on the thermohaline diffusivity in § 3.3. § 3.4 is de- 
voted to the case of stars in the mass range L5-2.2 Mq. Then 
in § 3.5 we shortly discuss the predictions for intermediate-mass 
stars. 



Figure |2] depicts the chemical structure of a 1.25 Mp star at 
the end of central hydrogen-burning (point Ai 25; top left panel 
for the present case without rotation-induced mixing). The most 
fragile elements lithium, beryllium, and boron, which burn at 
relatively low-temperatures and are preserved only in the most 
external stellar layers, are not shown here. On the pre-main se- 
quence, pristine deuterium is converted to ^He, while on the 
main sequence H-burning through the pp-chains builds up a ^ He 
peak at Mr/M* ^^0.65. Deeper inside the star the '^C peak re- 
sults from the competition between the '^C(p,7)'''N(j3 + v)'^C 
and '^C(p,7)'^N reactions. '^C and '^O are partially converted 
into '^N which abundance profile presents a double plateau. ON- 
cycling results in ^^O depletion and in the building up of a '^O 
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log (Teff) MJM, My'H. 



Fig. 1. Evolutionary track (from the pre-main sequence up to the 
tip of the RGB) of the 1.25 M^^ model computed with thermoha- 
line mixing only. The minimum and maximum luminosity of the 
bump are indicated (Lb.min and Lb.max respectively) as well as the 
luminosity th at which the thermohaline instability reaches the 
bottom of the convective envelope. Open symbols labelled Ai 25 
to El 25 correspond to evolution points for which some stellar 
properties are discussed in the text. The panel inserted on the 
right of the figure shows the evolution of the stellar luminosity 
around the bump as a function of time. Atb is the time spend by 
the star within the luminosity bump and equals to 3.9 x 10^ yrs 
in the present case. 



peak. In the very central regions, Na is produced through pro- 
ton capture by ^^Ne. 

When the star moves towards the RGB its convective enve- 
lope deepens and engulfes most of the regions that have been nu- 
clearly processed (IDUP). In Fig.|2]the maximum depth reached 
by the convective envelope is indicated by the vertical arrow. The 
so-called first dredge-up results in severe changes in the surface 
chemical properties of the star (see Fig.|6l) since surface material 
is diluted with matter enriched in ^He, '^C, and '"^N, but depleted 
in '^C and '^O. In the standard models the surface abundances 
are not predicted to change anymore after the end of the first 
dredge-up until the star reaches the early- AGB. However, as we 
shall see below, thermohaline mixing induces a second modifi- 
cation of the stellar chemical appearance on the upper end of the 
RGB. 

3.1 .2. Predictions up to tine RGB tip 

As discussed in § 1, the thermohaline instabiUty induced by the 
^He(^He,2p)^He reaction is able to set in on the RGB only after 
the star has reached the luminosity bump, when the HBS crosses 
the molecular weight barrier left behind by the convective enve- 
lope when it reached its maximum extent. In the case where ther- 
mohaline mixing is the only transport process considered within 
radiative regions (i.e., rotation-induced mixing being neglected), 
we find that for solar metallicity stars with initial mass lower or 



Fig. 2. Chemical structure at the turnoff of the 1 .25 star com- 
puted for different initial rotation velocities as indicated. The 
mass fractions are multipHed by 100 for ^He, ^^c, and ''^N, by 
2500 for "C, by 50, 900, and 5 x 10"* for '''O, "O, and '^O re- 
spectively, and by 1500 for ^^Na. The vertical arrows show, in all 
cases, the maximum depth reached by the convective envelope 
at its maximum extent during the first dredge-up. 

equal to 1.5 M,:;, the thermohaline instability soon extends be- 
tween the external wing of the HBS and the base of the con- 
vective envelope (see Table 1). This is shown for the 1.25 M© 
model in Fig. |3] where the thermohaline region is indicated in 
blue. For this model, the bump is located between Lb,min=37 Lq 
and Lb,max=45 Lq, and the thermohaline instability contacts the 
convective envelope when the stellar luminosity is Lc th ~ 94 Lq. 
As far as timescales are concerned, the 1.25 Mq model spends 
3.9x10^ years in the bump (i.e., in the luminosity interval be- 
tween Lb min and Lb,max> See Fig. [Til, and then reaches Lc.th after 
2.65.10^ years. 

Figure|4]depicts the temporal evolution along the RGB of the 
abundance profiles of H, ^He, of the carbon isotopic ratio, and of 
the thermohaline diffusion coefficient Dthc- The variable 8M is a 
relative mass coordinate allowing for a blow-up of the radiative 
region above the HBS. 

Namely, 

Mbce - Mubs 

It is equal to 1 at the base of the convective envelope Mbce 
and to at the base of the HBS Mhbs (which is defined as the 
depth where the hydrogen mass fraction equals 10^'") . On each 
graph the black solid lines correspond to the evolution point 
B 1 25 on the HRD when the thermohaline region is still quenched 
in a very tiny region, while the other curves correspond to lat- 
ter times (points C1.25, D1.25, and Ei 25) when the thermohaline 
instability has extended up to the base of the convective enve- 
lope. The maximum depth of the thermohaline unstable region 
corresponds to the layer where the ^He(^He,7)^Be reaction be- 
comes faster than ^He(^He,2p)^He (5M^0.2; see upper panel in 
Fig. |5]l. There the hydrogen profile shows a peak which maxi- 
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0.2 - 



t(yr) 

Fig. 3. Kippenhahn diagram for the 1 .25 Mq star computed with 
thermohahne mixing (no rotation-induced processes). Here we 
focus on the evolution phase around the RGB luminosity bump 
(located at a total stellar luminosity of ^ 45 as indicated by 
the left arrow). Hatched area is the convective envelope. Green 
dotted lines delimit the hydrogen-burning shell above the degen- 
erate helium core. The zone where the thermohaline instability 
develops is shown in hatched blue. The right arrow indicates the 
total stellar luminosity (Lc,th ^ 94 L,7,) at which the thermoha- 
line region extends up to the convective envelope and thus when 
the surface chemical composition starts changing. 



mum is located at the depth where the reaction ^He(^He,2p)^He 
is the fastest. As soon as the thermohaline instability sets in, 
fresh protons diffuse outwards, spreading out the molecular- 
weight inversion and enlarging the thermohaline region until it 
reaches the convective envelope. Simultaneously ^He diffuses 
from the convective envelope inwards, which fuels the thermo- 
haline instability. '^C and '^C diffuse respectively inwards and 
outwards, leading to a decrease of the surface carbon isotopic 
ratio. '^N also diffuses outwards. Among the oxygen isotopes, 
only '^O is affected, which leads to a slight increase of the 
16q/18q suj-face ratio while '^O/'^O does not change. Elements 
with higher atomic numbers such as ^■'Na, which burn or are 
produced at higher temperature than that of the maximum depth 
of the thermohaline region, are unaffected. On Fig|5] one can 
see also that in the external wing of the HBS below 5M ~ 0.6, 
the thermohaline diffusion timescale T(diff) is longer than both 
T(^Be + e^) and T(^LiH-p). As a consequence no fresh ^Li shows 
up at the stellar surface on the RGB for the value of the coef- 
ficient Ct chosen for the computations. Rather, the surface ^Li 
abundance decreases as this fragile element is drained from the 
convective envelope downwards. For the 1.25 M,3 non-rotating 
model the value of N(Lij3at the RGB tip is - (see Table 2 and 
Fig.[TO]in § 3.3). 

The evolution of the surface abundances from the zero age 
main sequence up to the RGB tip is shown on Fig|6](left pan- 
els for the present case without rotation-induced mixing). One 
sees clearly the signatures of the first dredge-up at Log(L/L©) ^ 
0.6 that leads to an increase of the surface abundances of ^He, 
'^C, '"^N, and '^O, and of the ^''O/'^O ratio, and to a decrease of 



N(Li)=log[n(Li)/n(H)]+12 




O.Z 0.4 0.6 0.8 10 0.2 0.4 0.6 0.8 1 
<5M 6U 

Fig. 4. Profiles of the abundances (in mass fraction) of H, ^He, 
of the carbon isotopic ratio, and of the thermohaline diffusion 
coefficient at various evolution points on the RGB (see Fig.[T]i: 
Slightly above the bump (point B1.25, black solid line), at the 
luminosity of the contact (point Ci 25, red dashed line), close 
from the RGB tip (point Di 25, magenta long dashed line), and 
at the RGB tip (point E1.25, blue dot-dashed line). The abscissa 
is the scaled mass coordinate 5M that allows a blow up of the 
region of interest (5M=0 at the base of the HBS and 5M=1 at 
the base of the convective envelope). 

'^C, '^O, of the '^c/i^Cand '''0/"0 ratios as the convective en- 
velope digs into the regions that have been nuclearly processed 
on the main sequence (see Fig|2]i. Then all surface abundances 
level off while the star ascends the RGB up to the bumf0. At 
log(L/L0)^2, thermohaline mixing does lead to a second change 
in the stellar surface composition. It induces in particular a sec- 
ond drop of the carbon isotopic ratio and efficiently destroys ^He 
and ^Li (see Fig. [TOl i while '"^N and '^O/'^O slightly increase. 
We note that the ^He surface abundance at the tip of the RGB 
is much reduced compared to its value after the first dredge-up, 
although it remains higher than the initial value the star is born 
with. This is due to the combination of several factors related 
to the thermohaline diffusion timescale on one hand, and to the 
compactness of the HBS. More specifically, the final ^He surface 
abundance can in principle decrease down to the value of ^He at 
the bottom of the thermohaline unstable region (i.e., at 5M'^0.3 
in the 1.25 M© case; see Fig. IH. However, the timescale for 
thermohaline diffusion in this model is such that the surface ^^He 
at the RGB tip saturates at a higher value. In the case of low- 
mass, low-metallicity stars presented in CZ07, the thermohaline 
unstable region is more compact and has a steeper temperature 
gradient, resulting in a more efficient decrease of the surface ^^He 
abundance than in the present case. The same reasoning applies 
to the surface abundance changes in '^C and ^^C so that the car- 
bon isotopic ratio at the tip of the RGB saturates here to a value 

^ In the standard case without thermohaline mixing, the abundances 
obtained at the end of the first dredge-up remain unchanged until the 
convective envelope deepens for the second time on the early- AGB. 
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log(L/La) 



Fig. 5. 1 .25 M;7, model with thermohaline mixing (no rotation) at 
the evolution point C1.25 - (top) Lifetime of ^He, ^Li, ^Be, and 
'^C, and convective and thermohaline diffusive timescales (in 
years). The transition between the convective envelope and the 
radiative region occurs at 5M=1. (bottom) Abundance profiles 
(in mass fraction) of relevant species in the same region. 



of the order of 10, while lower ratios closer from the equilibrium 
value are reached in low-mass, low-metallicity stars (see Fig. 3 
and4ofCZ07). 

3.1 .3. Predictions on tine AGB 

After helium ignition at the RGB tip, the total stellar luminosity 
drops down to the location of the clump (at Log(L/L0) 1 .7 for 
the 1.25 Mf, model discussed here). A very modest decrease of 
the surface ^Li and '■^C/'^C occurs when the star starts ascending 
the early-AGB due to a second deepening of the convective en- 
velope (see Table 2). Also, ^He (in mass fraction) decreases very 
slightly from 3.7.10^'* to 3.3.10^^. In other words, some ^He re- 
mains in the external convective layers when the star enters the 
thermal pulse phase on the AGB (TP-AGB) and is available to 
drive thermohaline mixing further as discussed below. 

This 1 .25 Mq model was computed until the end of the su- 
perwind phase (at that point the total stellar mass is O.567M0 
and the mass of the envelope is O.OI7M0). Meanwhile it has 
undergone five thermal pulses. In the present study we assume 
no convectively induced extra-mixing below the convective en- 
velope as usually required to induce the third dredge-up in 
low-mass TP-AGB models jHerwig et alJll997HHerwidl2000l: 



Herwig et al l2007HMowlavill999H Karakas et alj2002tlKarakasl 
2010h . Consequently we do not expect our models to mimic 
carbon-rich stars (but see § 3.4). Such a study is out of the 
scope of the present work and TP-AGB predictions based on 
the present computations but including also different overshoot 
prescriptions will be presented in a future paper. 

Let us note, however, that we find that with the present as- 
sumptions thermohaline mixing is active during each interpulse 
and modifies the surface abundances as can be seen on Fig|2l As 



Fig. 6. Evolution of the surface abundances of ^He and '^N (in 
mass fraction) and of the carbon and oxygen isotopic ratios as 
a function of stellar luminosity for the 1.25 Mq models com- 
puted without (left) and with rotation (initial rotation velocity of 
110 km.sec^', right). Thermohaline instablity is accounted for 
in both cases. In the left panels the black and red lines corre- 
spond to Ct = 10^ and 10^ respectively, while in the right pan- 
els Ct = 10^. Predictions are shown from the zero age main se- 
quence up to the RGB tip. In the rotating case the first dredge- 
up starts at lower luminosity and ends at log(L/L0)'--^I with 
slightly stronger abundance modifications than in the standard 
case. Also, thermohaline mixing sets in slightly earlier on the 
RGB at log(Lc.th/L(; )'^l -9 when rotation is accounted for. 



on the RGB, thermohaline mixing changes, but in a very modest 
way, the abundances of ^He, '^C, '^C, "^N, ^^O, and '^O, while 
heavier elements are unaffected. 

Also and in contrast with what happened on the RGB, the 
surface abundance of ^Li now increases (see Table [3]). This is 
due to the mixing efficiency that allows the transport of fresh 
^Be outwards to regions cool enough for ^Li to survive. In the 
present case the production of ^Li is significant, and the surface 
N(Li) reaches a value of ~0.9 (then the star has a total lumi- 
nosity Log(L/L0) between ^3.0 and 3.6). We thus confirm, but 
this time at solar metallicity, the finding by ,Stanclif fe (2010) that 
thermohaline mixing does increase the surface Li abundance in 
low-mass TP-AGB stars. Whether these Li-rich objects can si- 
multaneously be carbon-rich will require further investigation of 
the TP-AGB phase including parametric convective overshoot 
as mentioned before. As can be seen in Table [3] the total stellar 
yields for lithium remain negative. 

3.2. 1.25 M0, model with tiiermolialine and rotation-induced 
mixing 

Let us now discuss the case of a 1.25 M0 model that does take 
into account both thermohaline instability and rotation-induced 
mixing as described in § 2.2 and 2.3. 
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Time since 1" TP (yr) 



Fig. 7. Evolution of the surface abundances of ^He and ^Li {top, 
dashed and solid lines respectively), of '"^N and of the carbon 
isotopic ratio (middle, solid and dashed lines respectively), and 
of the helium-burning luminosity and total stellar mass (bot- 
tom, solid and dashed lines respectively) on the TP-AGB of the 
1.25 Mq model computed with thermohaline mixing only. The 
abscissa is the time since the first thermal pulse. 



Fig. 8. Evolution with time along the RGB of the effective tem- 
perature and of the luminosity of the 1.25 Mg models computed 
with thermohaline mixing only (black solid line) and with both 
thermohaline and rotation-induced mixing (initial rotation ve- 
locity of 1 10 km.s^' 

, red dotted line). 



3.2.1. Main sequence and subgiant branch 

As described above, the thermohaline instability induced by 
^He-burning sets in only on the RGB after the bump. Before a 
star reaches that phase, the modifications of its internal and sur- 
face chemical abundances are thus driven by rotation-induced 
mixing on the main sequence and convective dilution during 
the first dredge-up episode on the subgiant branch and early- 
RGB. The predictions of our rotating models up to that phase 
have been extensively tested in previous papers. They account 
nicely for the behaviour of lithium and beryllium at the sur- 



^ 4.1; Talon & Charbonnell 19981 201C 


; Charbonnel & Talon 


1999'; 'Palacios et al."2003'; 'Pasquini et al. 


120041 Smilianic et al. 


2009; Charbonnel & Lagarde 2010). 



Rotation-induced mixing has also an impact on the internal 
abundance profiles of heavier chemicals involved in hydrogen- 
burning at higher temperatures than the fragile Li and Be. This 
can be seen at the moment of the turnoff in Fig.|2]for the 1 .25 Mg 
star computed for different initial rotation velocities. 

In the rotating models, the abundance gradients are 
smoothed out compared to the standard case: ^He, '^C, '^N, 
and ^^O diffuse outwards, while ^^C and '^O diffuse inwards. 
However, rotation-induced mixing is not efficient enough to no- 
ticeably change the surface abundances of these elements while 
on the main sequence for the 1 .25 M© modefl although it sets 
the scene for abundance variations in latter evolution phases. 

In particular the surface abundance variations during the first 
dredge-up are slightly strenghtened when rotation-induced mix- 



The only exception is Li which is strongly depleted in the rotating 



ing is accounted for, as shown in Fig. |6] For example, more 
^He is brought into the stellar envelope, and the post dredge- 
up '^C/'^C and '^O/'^O ratios are lower than in the non-rotating 
case. We note, however, that neither '^O nor ^^Na are affected, 
and that no ^^Na enhancement is expected at the surface of such 
a low-mass star during the first dredge-up. 



3.2.2. Red giant brancli 

The modifications of the internal abundances due to rotation- 
induced mixing that we just discussed for the 1.25 M,;; model 
do induce slight variations of the overal stellar structure and 
of the evolutionary track. In the case of such a low-mass star, 
the impact on the effective temperature and luminosity along 
the RGB is relatively modest as can be seen in Fig. [8] We note 
that in the case of the computation without rotation (black full 
line) the drop in luminosity at the RGB bump (that is associ- 
ated with a slight increase in Teff at ^ 5.8 x 10 years) is clearly 
noticeable. In the rotating case the inflexion in luminosity at the 
bump is less pronounced and occurs at slightly lower luminosity. 
Consequently, the stellar luminosity at which the thermohaline 
instability reaches the convective envelope is also slightly lower 
in the rotating case (see Table 1). 

Figure |9] shows, at the evolution point Ci 25 for the 
1.25 Mq model computed with an initial rotation velocity of 
1 10 km.sec^', the diffusion coefficient associated to the thermo- 
haline instability, Dthc (Eq. 1 for Ct = 1000), and the total dif- 
fusion coefficient associated to rotation, Dioi, that characterizes 
the transport of chemicals throu gh the i nterac tion of meridional 
circulati on and shear turbule nce IZalinl ( Il992i see e.g. Eqs. 5, 7 
and 8 in lPalacios et al.ll2006h . Dthc is five to six orders of magni- 
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Fig. 9. Thermohaline diffusion coefficient Dthc (black solid line) 
and total rotation coefficient Djot (red dashed line) as a function 
of 5M for a 1 . 25 Mq star (initial rotation velocity of 1 1 km s ^ ' ) 
on the RGB at a total stellar luminosity of 105L© (i.e., at the evo- 
lution point Ci.25)- In the current model at this precise evolution 
point the size of the region between 5M=0 and 1 is 0.687696 
while the thermohaline instability extends over 0.66136 R©. 



tude higher than Drot, and this result is independent of the initial 
rotation velocity on the ZAMS. 

T his confirms the finding by ICantiello & Langeil (l2008l 
that in the relevant layers thermohaline mixing has much 
higher diffusion coefficients than rotational and magnetic in- 
stabilitiejZI . I t is al so perfe ctly consistent with t he results of 
iPalacios et al.l ( l2006l see also lChaname et ani2005h who studied 
the impact of rotation-induced mixing on the RGB for low-mass 
stars and showed that it cannot (with the present assumptions) 
account for the abundance anomalies observed in bright giants. 
These authors noted that assuming differential rotation (i.e., uni- 
form specific angular momentum) instead of solid body rotation 
(i.e., uniform ang ular velocity) in the con vective envelope along 
the RGB (see, i.e. jBrun & Palaciosl2009h does lead to higher ef- 
ficiency of the rotation-induced processes below the convective 
envelope. However, even in that case, Palacios and collaborators 
showed that the total transport coefficient associated to rotation 
does not rise above 10^ cm^ s^' in the outer HBS, which is still 
much lower than the thermohaline diffusion coefficient. 

Thermohaline mixing thus governs the surface abundance 
variations on the upper half of the RGB as already discussed 
by CZ07. The corresponding predictions for the 1.25 M© model 
computed with an initial rotation velocity of 1 10 km.s^' can be 
seen in Fig.|6](right panels). 



^ This conclusion on the magnetic diffusivity was obtained in the case 
of magnetic fields that are created in differentially rotating star and is 
not valid for magnetic stars that poss ess anomalous fossil fields s uch as 
the Ap star descendants discussed in lCharbonnel & ZahnI ( 12007 ah . 



3.2.3. AGB 

The 1.25 M(7, model with an initial rotation velocity of 110 km 
s^' was computed until the end of the superwind phase (total 
stellar mass and mass of the envelope being respectively equal 
to 0.566 Mf , and 0.014 Mg), and has undergone four thermal 
pulses. During the TP-AGB the behaviour of the surface abun- 
dances is similar to that discussed in § 3.1.2. In the present case, 
the maximum N(Li) value reached at the end of the TP-AGB is ~ 
0.8 (instead of 0.9 in the 1.25 M ^) non-rotating model discussed 
in § 3.1; see also Table|3]l. Again, predictions for carbon at that 
phase must be taken with caution since the impact of parametric 
convectively induced extra-mixing is not taken into account. 

3.3. Uncertainties on fhe tliermolialine diffusion coefficient 

3.3.1 . Size and shape of the thermohaline cells 

CZ07 performed computations for various values of the coeffi- 
cient Ct and discussed the uncertainties on the efficiency of the 
thermohaline mixing that are basically related to the size and 
shape of the thermohaline cells. Their prefered value for the as- 
pect ratio a=5 (see § 2.2) also used in the present computations 
corresponds to "fingers" rather than "blobs" whose shorter mix- 
ing length would translate into smaller value (by a factor of ~ 
50) for the coefficient Ct . Crude as it may be, this choice first ad- 
vocated bv lUMchl (Il972h is supported by laboratory experiments 
where the instability takes the form of salt fingers ( Krishnamurt^ 
l2003h . Also and contrary to the "blob assumption", the "fin- 
ger prescription" leads to a very good description of the surface 
abundances in low-metallicity stars as shown in CZ07, as well as 
in solar-metallicity stars as will be discussed in § 4. Unfortunatly 
an d as mentionne d by Eggleton et al.l (|2007|) the 3D simulation 
bv lEggleton et al.l (l2006[i did not have the resolution to give clues 
on the aspect ratio of the fingers. 

As a test we have however run a 1.25 M^^:, model without 
rotation-induced mixing but with Ct = 10^ instead of the value 
of 10^^ used in all the other computations presented in the present 
paper The predictions for the evolution of the surface abun- 
dances up to the RGB tip are shown in Fig|6] (red lines in left 
panels) and[TO] As expected, increasing the thermohaline diffu- 
sion coefficient by a factor of 10 leads to faster and substantially 
stronger processing of material on the RGB: At the RGB tip, the 
surface abundances of ^He and of ^Li are reduced respectively 
by a factor of 2 and by ~ 1.5 dex compared to the Ct = 10^ as- 
sumption. The impact on the carbon and oxygen isotopic ratios 
and on the surface abundance of ''^N is more moderate. 

Even in that case there remains enough ^He to drive ther- 
mohaline mixing when the star is on the TP-AGB. The Li pro- 
duction during that phase is higher than in the Ct = 10^ case, 
with the final surface abundance N(Li)=2 instead of 0.9. The fi- 
nal ^He abundance and carbon isotopic ratio are 1.8 x 10^^ (in 
mass fraction) and 8 respectively (instead of 3.15 x 10^^ and 
9.65). 

3.3.2. Atomic diffusion 

In the present computations we have not included the effects of 
atomic diffusion. In particular we do not consider radiative levi- 
tation that may lead to accumulation of heavy elements and thus 
to thermohaline instability in the outer layers of peculiar, slowly- 
rotating main sequence A-type stars (see Theado et al. 2003)- 
This simplification has no effect on our conclusions, since the 
thermohaline instability induced by iron accumulation affects 
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Fig. 10. Evolution of the surface abundance of ^Li in the 
1 .25 Mq star up to the RGB tip when considering thermohaline 
transport but no rotation-induced mixing. The black solid and 
red dashed curves correspond to computations performed with 
Ct = 10"* and 10^ respectively. 



Fig. 11. Same as Fig. |2]for the 2.0 M,7, star computed for dif- 
ferent initial rotation velocities, as indicated. The mass fractions 
are multiplied by 100 for ^He, '^c, and ''^N, by 2500 for '^C, by 
50, 1 100, and 5 x 10"^ for '*^0, "O, and ^^O respectively, and by 
1500 for ^^Na. The vertical arrows show, in all cases, the maxi- 
mum depth reached by the convective envelope at its maximum 
extent during the first dredge-up. 



only the very external regions of these atypical main sequence 
stars and has no direct impact on the nuclear burning occuring 
much deeper inside the star, nor on the RGB chemical properties. 

We do not con sider either the effect of atomic diffusion on 
the RGB, although lMichaud et all (1201 Oh pointed out that at that 
phase ^He gravitational settling may eventually lead to a larger 
/i-inversion than ^^He-burning on the outskirts of the HBS. In 
their computations, however, thermohaline mixing is not taken 
into account. Consequently the effects of concentration varia- 
tions on jj. they get from pure atomic diffusion are maximum 
compared to reality where thermohaline mixing (induced by 
^He-burning and eventually by '^He-settling) counteracts atomic 
diffusion. Michaud and collaborators have estimated that a value 
of Dthc of the order of 10^ cm^ s^' is able to substantially reduce 
(by a factor of 10) the small gradients of He caused by atomic 
diffusion on the RGB for a 0.95Me, Z=0.004 model. Given that 
this number is smaller than Dthc obtained in our RGB models 
(see Fig.|4]and|9l), we can safely assume that the effects of atomic 
diffusion must be wiped out by turbulence and that ^^He-burning 
is the dominating process inducing thermohaline instability be- 
tween the HBS and the convective envelope in RGB stars. We are 
aware that some "^He settling may remain even under the coun- 
teracting action of thermohaline mixing, although this should 
be confirmed by computations that are out of the scope of the 
present paper. One may note, however, that this should simply 
slightly re-inforce the /i-inversion induced by ^^He-burning (al- 
though to a much lower extent than in Michaud's computations), 
and thus strengthen the thermohaline transport compared to the 
present models. 



3.4. Low-mass stars more massive than ^7.5 Mq 

Table 1 gives the luminosity of the bump as well as the lu- 
minosity Lc^th at which the thermohaline instability "contacts" 
the base of the convective envelope for all the low-mass stellar 
models computed both without and with rotation for the present 
study. As we have just seen, for RGB stars less massive than 
1 . 5 M,7, , thermohaline mixing starts changing the surface abun- 
dances soon after the bump. However, for RGB stars with ini- 
tial mass higher than 1 .5 Mfr) computed without rotation-induced 
mixing, the thermohaline instability is long quenched into a very 
thin region, and is able to connect the external wing of the HBS 
with the convective envelope only when the star reaches already 
very high luminosit y, close from the RGB tip. This is consistent 
with the finding bv ICantiello & Langed (l2008i 1201 Oh of an up- 
per mass limit for efficient thermohaline mixing in non-rotating 
low-mass RGB stars. However, as we shall see below, this con- 
clusion does not hold anymore when one considers the impact 
of stellar rotation. 

Let us now discuss indeed the case of a 2.0 M,;j star which in- 
terior chemical structure at the turnoff is shown in FiglTT]for dif- 
ferent initial velocities. As for the 1.25 Mq models discussed be- 
fore, rotation-induced mixing smoothes the abundance profiles 
inside the star and in the present case it already leads to varia- 
tions of the surface abundances on the main sequence that are 
stronger for higher initial rotation velocities (this can be seen by 
looking at the abundance v alues at Mr/M^=l in the v arious pan- 
els of Figinj (see also e.g. lMevnet & MaedeJl2002h '). Note also 
that the maximum value of ^He at the peak is much lower than 
in the 1.25 Mq, because on the main sequence the 2.0 Mq, burns 
hydrogen mainly through the CNO cycle rather than through the 
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Fig. 12. Evolution of the surface abundances of ^He, ^Li, '^N, 
and of the carbon and oxygen isotopic ratios for the 2.0 M© 
star from the zero age main sequence up to the RGB tip. {Left) 
Model including thermohaline mixing only. (Right) Models in- 
cluding thermohaline and rotation-induced mixing for different 
initial rotation velocities (110 and 250 km s^', solid black and 
dashed blue lines respectively). In all cases Ct — 10^. 



Fig. 13. Chemical structure at the turnoff of the 4.0 M,7, star com- 
puted without rotation and with rotation for an initial rotation ve- 
locity of 300 km s^' (top and bottom panels respectively). The 
mass fractions are multiplied by 600 for ^He, by 100 for 
and i^N, by 2500 for '^C, by 50, 5000, and 6x 10"^ for '^O, '^O, 
and '^O respectively, and by 1500 for ^^Na. The vertical arrows 
show, in both cases, the maximum depth reached by the convec- 
tive envelope at its maximum extent during the first dredge-up. 



pp-chains. The more massive star thus dredges less ^He on the 
subgiant branch as can be seen by comparing Fig.l6landfT2] 

FigurdT2l shows the evolution of the surface abundances of 
^He, ^Li, i^N, and of the ^^C/^^C, ^^0/"0, and ^^O/^^O ratios 
as a function of luminosity up to the RGB tip in 2.0 Mq models 
computed without or with rotation (left and right panels respec- 
tively). At this stellar mass some ^^Na produced during the main 
sequence is dredged-up (not shown in Fi dT2l but see Table 2). 
All these quantities at the end of the first dredge-up are slightly 
affected by rotation-induced mixing that changed the abundance 
profiles while the star was on the main sequence (e.g., the post 
dredge-up value of the carbon isotopic ratio is lower when faster 
rotation is accounted for). Furthermore, we note that the changes 
in surface abundances due to thermohaline mixing start at much 
lower luminosity (closer from the bump) on the RGB than in the 
non-rotating case (see also Table 1). Overall, the total abundance 
variations at the tip of the RGB are stronger for higher initial ro- 
tation velocity. 

Th us and contrary to the conc lusion by ICantiello et al.l 
(12001) and lCantiello & Langeii(l2010l) . we find no mass limit (in 
the case of low-mass stars that ignite He in a degenerate core) for 
thermohaline mixing to change the surface abundances on the 
RGB. However, the global efficiency of this process increases 
when one considers less massive stars at a given metallicity, or 
more metal-poor stars at a given stellar mass. As discussed pre- 
viously, this results from the combination of several factors like 
the thermohaline diffusion timescale compared to the secular 
timescale, the compactness of the HBS and of the thermohaline 
unstable region, and the amount of ^He available to power the 
thermohaline instability. 



In this mass range thermohaline mixing also leads to Li pro- 
duction during the TP-AGB (see Table |3j. In the 2.0 Mj model 
computed with an initial rotation velocity of 1 10 km s^', N(Li) 
increases from a value of 0.1 on the early AGB up to 1.5 at 
the end of the TP-AGB (the envelope mass is then ~ 0.14 Mq, 
and the star has undergone 1 1 thermal pulses). Meanwhile the 
^He surface abundance has slightly decreased from 2.74 x 10^^ 
to 2.67 X lO-'*. 

We wish to emphasize an interesting result obtained for the 
2.0 Mg, model that was computed up to the AGB tip with ther- 
mohaline mixing but without rotation. This model did undergo 
11 thermal pulses in total. After the 9th thermal pulse third 
dredge-up occuiTed, that slightly increased the '^C surface abun- 
dance as well as the carbon isotopic ratio. During the following 
interpulse phase this ratio was slightly lowered under the effect 
of thermohaline mixing. Quantitatively, the carbon isotopic ra- 
tio increased from 19.6 to 21.33 between the end of the second 
dredge-up and the AGB tip (see also FiglTS] and discussion in 
§ 4.2). This could indicate that thermohaline mixing does favour 
the occurrence of third dredge-up. This will be investigated fur- 
ther in a future paper. 

3.5. Theoretical predictions for intermediate-mass stars 

By definition, intermediate-mass stars are objects that ignite cen- 
tral helium-burning in a non-degenerate core at relatively low 
luminosity on the RGB, well before the HBS reaches the mean 
molecular weight discontinuity caused by the first dredge-up. In 
other words, these objects do not go through a bump on their 
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short ascend of the RGB, and thus do not undergo thermohaline 
mixing at that phase. 

We should note, however, that as in the previous cases, 
rotation-induced mi xing can not be neglected from the whole 
picture. We refer to lEggenberger et aP (1201 (j) for a discussion 
of the global effects of rotation on the evolution and astero- 
sismic pr operties of intermediate-mass red giants (see also e.g. 
iMevnet & Maeder _2QQ2). Its impact on the chemical structure of 
a 4.0 M© star at turnoff can be seen in Fig. [13] Note that in this 
mass range the base of the convective envelope reaches the ^^Na 
plateau during the first dredge-up, leading to an increase of the 
surface abundance of this element both in the non-rotating and 
rotating cases. 

Overall rotation-induced mixing leads to stronger modifica- 
tions of the stellar chemical properties when the star becomes a 
giant as its convective envelope dredges-up nuclearly processed 
material. At the end of the dredge-up for the 4 M© models, the 
surface ■'He abundance is 1.1 x 10^"* and 8.4 x 10^^ in the non- 
rotating and rotating (initial rotation velocity of 300 km s^') 
models respectively, while the carbon isotopic ratio is 20.5 or 14 
respectively, and N(Li) is 1.3 or -1.8. 

We computed the first 1 1 thermal pulses for the 4 M© models 
without and with rotation-induced mixing, including thermoha- 
line mixing in both cases. A strong Li increase at the stellar sur- 
face is obtained during the first thermal pulses, and then lithium 
production levels off at a val ue of the order of N(Li)=2.2. This 
agrees with IStancliffd (1201 Ol) predictions. Again, a detailed ex- 
ploration of the TP-AGB phase for intermediate-mass stars is 
postponed to a further paper. 

4. Comparison with observations 

We now test the theoretical predictions of our models with re- 
spect to observations of relevant chemical elements in stars at 
different evolution stages. In Table |2] we give the surface car- 
bon isotopic ratio as well as surface abundances of ^Li, ^Be, 
[C/Fe], [N/Fe], and [Na/Fe], at the end of the first dredge-up, at 
the RGB tip, and at the end of the second dredge-up, for each of 
the models we have computed. In this table and in the following 
figures all the predictions correspond to models computed with 
a value of Ct ~ 10^ for thermohaline mixing (without or with 
rotation-induced mixing). As underlined in § 3, we did not in- 
clude parametric convectively induced extra-mixing during the 
TP-AGB phase so that our models are not expected to undergo 
third dredge-up and to mimic carbon-rich stars. 

4. 1 . Lithium beliaviour 
4.1.1. Lithium destruction 

As already mentioned the predictions of the present rotating 
models have been successfully compared to Li and Be obser- 
vations along the whole evolutionary sequence of the Galactic 
open c luster IC 4651 (turnoff mass 1.8 M©) bv lSmiljanic et al.l 
(1201 Ol see their Fig. 10 to 14). They account very nicely for all 
the Li and Be features observed in this cluster. 

Here we use as additional constraints Li observations that 
we performed for a large sample of field red giant stars (sub- 
giant, RGB, and early-AGB stars) with metallicities around so- 
lar. All sample stars have Hipparcos parallaxes so that their 
mass and evolutionary status could be relatively well determined 
(Charbonnel et al. in preparation). In FiglT4] they are distin- 
guished with respect to their mass (less or more massive than 
2 M© in the left and right panels respectively). 
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Fig. 14. Lithium data for field evolved stars from the sample by 
Charbonnel et al. (in prep., see the text) that are segregated ac- 
cording to their mass (left and right panels include respectively 
sample stars with masses lower and higher than 2 M,;,; Li detec- 
tions and upper limits are shown as circles and triangles respec- 
tively). Theoretical lithium evolution is shown from the ZAMS 
up to the end of the early-AGB. Various lines coiTespond to pre- 
dictions for stellar models of different masses computed without 
or with rotation as indicated, and with thermohaline mixing in 
all cases (with Ct = 10^). 



Let us consider first the stars with initial masses lower than 
2 M©,, whose Li properties are compared with predictions for 
the 1.5 and 2 M,;) models (left panel of FigfT4li. The theoretical 
Li behaviour is relatively straightforward. On the main sequence 
and on the early-RGB, rotation-induced mixing leads to stronger 
Li depletion than in the standard case (compare e.g. the red curve 
with the black one); 

in this mass range indeed standard models predict no Li de- 
pletion on the main sequence and a N(Li) of the order of 1.5 at 
the end of the first dredge-up, which is at odds with the data. 
After the end of the first dredge-up (Teff ^ 4800 K), the the- 
oretical Li abundance remains temporarily constant as the con- 
vective envelope withdraws in mass. When thermohaline mixing 
becomes efficient (Teff ^ 4200 K), the theoretical Li abundance 
drops again in drastic manner (while it would stay constant in 
the standard case). After the star has reached the RGB tip its ef- 
fective temperature increases (up to ~ 4800 K) as it settles on 
the clump, before decreasing again when the star starts climb- 
ing the early-AGB. The second dredge-up that occurs then leads 
to a final decrease of N(Li). On this graph we do not plot the 
Li increase that is predicted to occur during the TP-AGB phase 
at a Teff of ~ 3200 K due to thermohaline mixing, and which 
is discussed in § 4.1.2. As can be seen in FiglT4] the present 
predictions are in perfect agreement with the data all along the 
evolutionary sequence and explain very well the upper limits ob- 
served for the brightest sample giant stars. The observed Li dis- 
persion at a given effective temperature reflects dispersion in the 
initial rotation velocity and in the initial stellar mass (see also 
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1.25Ma Vj,„s=Okm.s-' 
1.25Ma V„„3=110km,s-' 
2.OM5, V„„3=0km.s-' 
2.0M. V,.„,= 110km.s-' 




Table 3. Surface lithium abundance after the second dredge-up 
and at the end of the TP-AGB phase, and total lithium yield. 
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Teff 
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Fig. 15. Lithium observations in lUttenthaler & LebzelteJ (1201 Ol) 
sample of oxygen-rich variables belonging to the Galactic disk 
(circles and triangles are for abundance determinations and up- 
per limits respectively) as a function of effective temperature 
and bolometric magnitude. Typical error bars are indicated. The 
star with the highest Li abundance is W441 Cyg (see text). 
Theoretical lithium evolution is shown from the early-AGB up 
to the end of the TP-AGB. Various lines correspond to predic- 
tions for stellar models of different masses computed without or 
with rotation as indicated, and with thermohaline mixing in all 
cases (with Ct = 10^). 



Charbonnel & Talonlll999l: iPalacios et al.ll2003l: ISmiljanic et al.l 
2010h . 



The case of the more massive stars, whose Li observational 
behaviour is compared to predictions for the 2.5 and 2.7 M© 
models (right panel of FigO is even more simple. In these 
objects indeed no thermohaline mixing occurs on the too short 
RGB, and rotation-induced mixing alone explains very well the 
data. 



4.1.2. Lithium production 

In all the models that we have computed along the TP-AGB, non 
negligible fresh lithium production is obtained, although the to- 
tal stellar yields remain negative. Table |3] gives the final N(Li) 
value and the total Li yield for those models. In all cases ther- 
mohaline transport is responsible for this strong Li enrichment. 

The evolution of the surface abundance along the TP-AGB 
is shown as a function of effective temperature and bolometric 
magnitude in FiglT5]for the 1.25 and 2.0 M© models computed 
without and with rotation, and with thermohaline mixing in both 
cases (with Ct = 10^). Li production starts at slightly higher ef- 
fective temperature and luminosity for the more massive star. 

Predictions are compared with lithium values in the sam- 
ple of low-mass oxygen-ric h AGB vari ables belonging to 
the Galactic disk studied by lUttenthaler & Lebzelter (2010). 
Theoretical Li production sets in above a lower luminosity limit 
which agrees with the observational Mbol threshold. Models are 
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Table 4. References for the abundance studies in Galact ic open 
cluster s used i n the compariso ns with model predictions 
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ISmiiianic et alJ (l2009l) : S09; lYong et al.l (12001: Y05. The red 



turnoff masses given in the second column were estimated using 
the WEBDA database and Geneva isochrones (see text). Listed 
in the last column are the symbol colours used in Fig. 17 and 
Fig. 19 to 21. 
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found to fit very nicely the lithium behaviour. However, they can 
not account for the very high Li abundance (3.1 to 4.6 depend- 
ing on the model atmosphere) of the star V441 Cyg, which may 
rather be an intermediate- mass AGB star undergoing hot bottom 
burning (see discussion in lUttenthaler & Lebzelteijr2010l) . 

4.2. Carbon isotopic ratio 

As discussed in the introduction, the behaviour of the carbon 
isotopic ratio is the best indicator of non-standard transport pro- 
cesses in evolved low-mass stars. This quantity has been de- 
termined in a large number of stars in Galactic open clusters. 
The references for the studies we consider here are listed in 
Table 4 where we also give the cluster red turnoff masses (see 



14 



C.Charbonnel & N.Lagarde: Thermohaline instability and rotation-induced mixing 




30 



1 2 
Log(L/Lj 

Fig. 16. Evolution of the surface '^C/'^^C value as a function 
of stellar luminosity for the 1.25M,;, models including thermo- 
haline instability and rotation-induced mixing (for initial rota- 
tion velocities of 50, 80, and 110 km.s^' shown as solid red, 
dashed green, and dotted blue lines respectively). The non ro- 
tating case is also shown (black solid line). Observations along 
the evolutionary seque nce of the open cluster M67 are from 
lGilrov&Br6wnl ( ll991h . The triangle is for a subgiant star for 
which only a lower value could be obtained, while black squares 
and red circles correspond respectively to RGB and clump stars. 



iMevnet etal1ll993h derived from the WEBDA d atabase when 
using the Geneva isochrones (ISchaller et al.lll992l) . 

Figure [16] displays observations of the '^C/'^C ratio in sub- 
giant, RGB, and clump stars of the open cluster M67 (t urnoff 
mass ^ 1.2 M© according to iGilrov & Brownl (Il99lh. and 
1.5 M p, according to WEBDA, see above) by Gilrov & Brownl 
The data are compared to the predictions of our 1.25M0 
models computed for three different initial velocities and includ- 
ing thermohaline mixing (the non-rotating case is also shown). 
We see that the theoretical and observational behaviours are in 
complete agreement all along the evolutionary sequence. While 
the dispersion for stars that have not yet reached the RGB bump 
(i.e., with Log(L/L)0 between ^ 0.7 and 1.8) reflects only the 
dispersion in the initial rotation velocity, explaining the data for 
more evolved stars requires the occurrence of thermohaline mix- 
ing as predicted by the models. 

We now compare the predictions of our models over the 
1.0-4.0 M0 range with carbon isotopic ratios in open clus- 
ters of different turnoff rnasses . The data sho wn in Fig. [iTl 
are from iGilrov & Brownl (Il99ll) for M67, from iGilrovl (1 19891) 
for her open clusters with turnoff mass es below 1.7 Mp, (i.e. , 
NGC 752), from the more recent study bv lSmilianic et all (l2009h 
for nine pen clusters with turnoff masses above 1 .7 Mq (IC 27 14, 
IC 4756, NGC 2360, NGC 2447, NGC 35 32, NGC 5822 
NGC 6134, NGC 6281, NGC 6633), and from lMikolaitis et alJ 
(I2OIOI) for NGC 6134. Individual stars are attributed the red 
turnoff mass of their host cluster determined as described above 
(see Table 4). Indications on their evolutionary status, when 
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Fig. 17. Observa tions of '^C/'^C in ev olved stars of Galactic 
open c l usters by Smilianic et al.l ( |2009[ open symbol s), Gilrovl 
(119891) . iGiirov & Brownl ( 1199 lb . and iMikolaitis et all ( l2010l) as 
a function of the turnoff mass of the corresponding host cluster 
that can be identified thanks to the colours of the symbols (see 
text and Table |4]i. Squares, circles, and asteriscs are for RGB, 
clump, an d early- AGB s tars respectively, while diamonds are for 
stars from lGilrovl (Il989l) sample with doubtful evolutionary sta- 
tus; triangles are for lower limits. A typical error bar is indicated. 
Theoretical predictions are shown at the tip of the RGB and after 
completion of the second dredge-up (black and blue lines respec- 
tively). Standard models (no thermohaline nor rotation-induced 
mixing) are shown as dotted lines, models with thermohaline 
mixing only (Vzams=0) as solid lines, and models with thermo- 
haline and rotation-induced mixing for different initial rotation 
velocities as indicated as long-dashed, dot-dashed, and dashed 
lines. 



available, are given in the plot ( squares, circles, and asteriscs are 
for RGB, clump, early-AGB stars respectively, while diamonds 
are for stars with uncertain evolutionary status). Model predic- 
tions are shown both at the tip of the RGB and at the end of the 
second dredge-upQ (black and blue lines respectively) for dif- 
ferent assumptions. Dotted lines correspond to standard models 
computed without thermohaline mixing nor rotation; those ac- 
count only for the upper envelope of the data. Solid lines corre- 
spond to models computed with thermohaline mixing only. They 
reproduce very well the '^C/'^C behaviour for stars with initial 
masses lower than ~1.7M0. In this mass range rotation-induced 
mixing leads only to slightly lower values as shown by the long- 
dashed lines for an initial velocity of 1 10 km s^'. Note that the 
squares at Mtumoff= 1 -5 M© correspond to the M67 stars that have 
not yet reached the RGB bump (see FigfTSll. which explains why 
they lie between the standard and the thermohaline curves. For 
stars with masses between ^-^1.7 and 2.2M0 both thermohaline 



^ Note that with the present computations without parametric con- 
vective overshoot during the thermal pulses, the carbon isotopic ratio is 
only very slightly modified during the TP-AGB phase compared to its 
value at the end of the second dredge-up. 



C.Charbonnel & N.Lagarde: Thermohaline instability and rotation-induced mixing 



15 




2 3 
Progenitor Mass (Mg) 

Fig. 18. Observations of ^^ C/'^C in planetar y nebulae as a func- 
tion of the progenitor mass ( iPalla et al.ll2000h . Squares and trian- 
gles are respectively for actual determinations and lower limits 
of the carbon isotopic ratio. Black and blue theoretical lines have 
the same meaning as in Fig[T7l while red curves show the model 
predictions at the AGB tip for the low-mass models that were 
computed up to that phase (solid and dashed red curves are re- 
spectively for the models with thermohaline mixing only and for 
the models with thermohaline and rotation-induced mixing). 



and rotation-induced mixing are required to fit the data. For more 
massive stars thermohaline mixing plays no role but the observa- 
tional uncertainties allow the data to be well accounted for when 
rotation-induced mixing is taken into account (dot-dashed and 
dashed lines for initial velocities of 250 and 300 km s^' respec- 
tively). 

Finally we compare in Fig|T8]our predictions for the carbon 
isotopic ratio to its determination in a sample of planetary nebu- 
lae obtain ed by means of m illimeter wave observations of '^CO 
and i^co (iPalla et al.ll20"00l) . The abcissa is the progenitor mass 
derived by Palla and collaborators; this quantity is highly un- 
certain as shown by the error bars. Since the formation of the 
planetary nebula occurs at the AGB tip, the data should be com- 
pared to the model predictions at the end of the superwind phase 
as plotted in red. As explained in § 3 only some of our low- 
mass models (with initial mass < 2 Mq) were computed up to 
that phase. For these objects thermohaline mixing was found to 
slightly lower the carbon isotopic ratio during the thermal pulse 
phase, except in the 2 model with thermohaline mixing and 
no rotation that has undergone third dredge-up from the 9th pulse 
on (see § 3.4); in that case the carbon isotopic ratio at the AGB 
tip is slightly higher than at the end of the second dredge-up. 
Overall the comparison between the models and the data turns 
out to be quite satisfactory. Two of the planetary nebulae with 
low-mass progenitors (namely NGC 6781 and M 1-17) actually 
exhibit relatively high carbon isotopic ratio. Although the uncer- 
tainty on the initial stellar mass of these objects allows the data to 
be well accounted for by the thermohaline models (both with and 
without rotation), it could also be that in this couple of stars ther- 
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Fig. 19. [N/C] ratio as a function of turn-off mass of the host 
cluster for the Galactic open c l uster s ample by ISmilianic et al.l 
( I2009h and for M67 bv lBrownl ( Il987h . Symbols and lines have 
the same meaning as in FiglTTl 



mohaline mixing was inhibited by strong fossil magnetic fields 
as suggested by ICharbonnel &Zahn (2007 a). In this context it 
would be extremely valuable to look for magnetic fields in these 
two possible "thermohaline deviant stars". On the other hand, 
computations are now needed to estimate the combined effect of 
third dredge-up, hot bottom burning, and thermohaline mixing 
during the TP- AGB phase for stars more massive than 2 M©. 
This is work is in progress. 

4.3. Nitrogen, sodium, and oxygen isotopes 

[N/C] data for the open cluster sample listed in Table 4 is shown 
as a function of the cluster red turnoff mass and as a function of 
'^C/'^C in Fig[T9land[20lrespectivelv. 

During the first dredge-up, the convective envelope of in- 
termediate mass-stars reaches the regions where '^O has been 
partially converted into '^N while the star was on the main se- 
quence, while in the case of low-mass stars it reaches only the 
first '^N step due to '^C-burning (compare the position of the 
vertical arrows in Fig.[2][TTl andfTSTl. Models thus predict an in- 
crease of the post dredge-up [N/C] value with initial stellar mass, 
in agreement with the observed behaviour. We note, however, 
that the lower envelope of the observational data lies slightly 
below the standard predictions, which might indicate that the 
models overestimate the first dredge-up. On the other hand, the 
corresponding offset may well be related to observational un- 
certainties. Over the whole mass range thermohaline mixing on 
one hand, and rotation-induced mixing on the other hand, lead 
to additional transport of CNO-cycled material, and thus further 
increase the [N/C] ratio with respect to the standard predictions. 
Given the observational error bars, one can conclude from FigfT9l 
and|20]that the models account nicely for the observational con- 
straints on C and N. 
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Fig. 20. '^c/i^C as a function of [N/C]. See Table 4 for refer- 
ences. Symbols are the same as in FigfTTl 
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Fig. 21. [Na/Fe] ratio as a function of turn-off mass for the sam- 
ple of Galactic open clusters listed in Table 4. Symbols and lines 
are the same as in Fig [17] The green diamon ds are the mean 
[Na/Fe] values given bv iJacobson et all (l2007l) for IC 4756 and 
NGC 6939. 



Fig. 22. "'O/'^O versu s "'0/"^0. Obser vations are from 
iHarris & Lambert! d 19841) . [Harris et all ( Il988l) for stars in the 1- 
3M range. The line symbols are the same as in Fig(T7l The 
initial values assumed for '^O/'^O and '^O/^^O are 1490 and 
445 respectively. 



position of our modelfl Both the predictions and observations 
show a positive correlation between [Na/Fe] values and stel- 
lar mass. Rotation-induced mixing leads to an increase of the 
amount of Na processed to the surface and allows an explana- 
tion for the observed dispersion. There is however an offset of 
about 0. 1 dex b etween the data and th e predictions as was al- 
ready noticed bv lSmihanic et al.l ( |2009|) who c ompared their ob - 
servations with standard model predictions bv iMowlavil d 1 9991) . 
As a matter of fact, very different observational Na abundances 
for giant stars have been reported in the literature. As can be 
seen in Fi g. 1211 some studies pr esent [Na/Fe] values as high as 
+0.6 dex (iJacobson et al. Il2007h. some only a mild overabun- 
dance of -hO.2 dex (iHamdani et al. I I2OOOI see Fig HB, and other 
solar values dSestito et al.l 2007 ). We refer to ISmihanic et al.l 
(I2009h for a discussion on the possible causes of these discrep- 
ancies. The present predictions for the stars more massive than 
^ 2 M© are actually in better a greement with the mild overabun- 
dance of [Na/Fe] measured bv lHamdani et al.l (l2000l) . 

Finally we show in Fig.felthe '^O/'^O v s ^^O /^^O for the 
G and K giants by [Harris & LambertI (11984 and [Harris et all 
([1988[) . Included in the figure are our predictions. As discussed 
in § 3, thermohaline mixing affects only slightly the '^O/'^O ra- 
tio, and leaves '^O/'^O unaffected; on the other hand, rotation- 
induced mixing lowers the '^O/'^O ratio, and helps account for 
the lowest "'0/"*0 values. Given the large observational uncer- 



In figure |2T| we plot the [Na/Fe] ratio for the open cluster 
sample listed in Table 4 as a function of cluster turnoff mass. 
Note that t he observational data were reported to the solar Na 
value from [Asplund et aP ([2005[) we assume in the initial com- 



^ ISmilianic etaP ( i2009h and iMikolaitis et"an l l2010h adopted 
solar abundances recommende d b y Grevesse et al. (1993^ (A(Na) 
= 6.33, A(F e)=7.50 ), whi le [jacob son et al. ( 2007, 2008, 2003) 
and Hamdani et al.l |200(]h used solar abundances given by 
[Anders & Grevessd ( il989h (•ArNa)=6.33'). In Fig. [21] the observa- 
tional data are reported to the sola r abondances values we use in our 
computations dAsplund et al.ll2005L A(Na)=6.20, and A(Fe) = 7.45). 
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tainties, the predictions are reasonably consistent with the O iso- 
topic ratios measured in RGB stars. 

5. Conclusions 

In the present paper we have investigated the effects of the ther- 
mohaline instability induced by ^He-burning that sets in above 
the RGB bump and of rotation-induced mixing on the evolu- 
tion and chemical properties of low- and intermediate-mass stars 
(1 to 4 Mq) at solar metallicity. All the stellar models were 
computed up to the end of the second dredge-up on the early- 
AGB, and some of them up to the end of the TP-AGB phase. 
Predictions are compared to data for lithium, '^C/'^C, [N/C], 
[Na/Fe], ^^'O/'^O, and i^O/^^O in giant stars with well-defined 
masses and evolutionary status on the RGB, clump, early-AGB, 
and planetary nebulae phases. 

We find that the theoretical and observational behaviours for 
these species are in very good agreement over the whole scruti- 
nized mass range. Thermohaline mixing is confirmed to be the 
main physical process governing the surface abundances of ^He, 
^Li, C, and N for stars more evolved than the RGB bump in all 
the models with initial masses below 2.2 Mq, although its ef- 
ficiency is increasing with decreasing initial stellar mass. In all 
cases ^He decreases by a large fraction in the stellar yields com- 
pared to the standard models, although we find that low-mass 
stars remain net producers of ^He (^He yields for stellar models 
over a broad range in both mass and metallicity will be published 
in a future paper). It is also found that thermohaline mixing leads 
to lithium production on the TP-AGB phase, as first shown by 
IStanchffd midt) in the case of low-metallicity stars. However, 
the Li yields remain negative, and these stars are not expected 
to contribute to Galactic Li enrichment. In one 2.0 Mq model 
computed up to the AGB tip thermohaline mixing was found to 
help initiating the occurrence of the third dredge-up during the 
TP-AGB phase. 

On the other hand, rotation-induced mixing modifies the in- 
ternal chemical structure of main sequence stars, although its 
signatures are revealed only later in the evolution when the 
first dredge-up occurs. It favours the occurrence of thermoha- 
line mixing in RGB stars in the mass range between '^1.5 and 
2.2 Mf ,. It accounts for the observed dispersion of abundances in 
stars of similar mass and evolutionary status, and is necessary to 
explain the features of CN-processed material in intermediate- 
mass stars. 

These results were obtained using the prescription for the 
turbulent dif fusivity related to the thermohaline instability advo- 
cated first bvlUlrichl(ll972h that is supported by laboratory exper- 
iments Krishnamurtil The same prescription was shown 
bv lCharbonnel & ZahnI (l2007bh to nicely account for the photo- 
spheric composition of low-mass, low-metallicity giant stars. 
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Appendix A: Nuclear reaction rates 

All reactions for hydro gen burning are coin puted with nominal 
NACRE reaction rates jAngulo et alJI 19991) . with the exception 
of: 

'^C(p,y) ^^N (IWiescher et al.l Il990l); '^C(p,n)'^N 
(lKoehler&0'brienlll989h : ^'^C{p,ay^B (ICaughlan & Fowled 



Il988h : '"^NCp, 7)1^0 (iMukhamedzhanov et alJ l2003h : 
2iNe(p,7)22Na, 22Na(p,7pNa 2^Na(p,a)2'' Ne, and 
3la(p,7)24Mg (llliadis et al.ll200ll) : 22Ne(p,7)23Na (iHale et all 
I2002h . 
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